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We consider the evolution of a scalar field propagating in Schwarzschild-de Sitter spacetime. The 
field is non-minimally coupled to curvature through a coupling constant ^. The spacetime has two 
distinct time scales, te ~ r^jc and tc — rc/c, where re is the radius of the black-hole horizon, re 
the radius of the cosmological horizon, and c the speed of light. When re ^ r^, the field's time 
evolution can be separated into three epochs. At times t <^tc, the field behaves as if it were in pure 
Schwarzschild spacetime; the structure of spacetime far from the black hole has no infiuence on the 
evolution. In this early epoch, the field's initial outburst is followed by quasi-normal oscillations, 
and then by an inverse power- law decay. At times t ^ tc, the power- law behavior gives way to a 



faster, exponential decay. In this intermediate epoch, the conditions at radii 



> 



and 



< 



both play an important role. Finally, at times t ^ tc, the field behaves as if it were in pure de 
Sitter spacetime; the structure of spacetime near the black hole no longer influences the evolution 
in a significant way. In this late epoch, the field's behavior depends on the value of the curvature- 
coupling constant ^. If ^ is less than a critical value ^c ~ 3/16, the field decays exponentially, with 
a decay constant that increases with increasing ^. If ^ > £,c, the field oscillates with a frequency 
that increases with increasing ^; the amplitude of the field still decays exponentially, but the decay 
constant is independent of ^. We establish these properties using a combination of numerical and 
analytical methods. 

Pacs numbers: 04.30.Nk, 04.40.-b, 04.70.Bw 



I. INTRODUCTION AND SUMMARY 

The dynamics of radiative fields in black-hole space- 
times has been investigated since the early 1970's in an 
effort to understand how asymmetric gravitational col- 
lapse proceeds, in the absence of rotation, to form a per- 
fectly spherical black hole. The first work of this type was 
carried out by de la Cruz, Chase, and Israel [Q, who fo- 
cused on electromagnetic and linear gravitational pertur- 
bations of a spherical collapse. They showed numerically 
that after an outburst of radiation at the onset of col- 
lapse, the radiation decays to zero in a time of the order 
of the light-crossing time of the resulting event horizon. 
Subsequently, Price f^ gave an analytical description of 
the radiative decay, and showed that scalar, electromag- 
netic, and gravitational radiation all go to zero as an 
inverse power of time, the power index depending on the 
multipole order of the perturbation. This late-time de- 
cay of the radiation is often referred to as its "power-law 
tail" . 

The early investigations were extended and improved 
by other authors. Bicak g generalized Price's analysis 
and calculated the radiative decay of an electrically neu- 
tral scalar field interacting with a charged black hole. He 



also found an inverse power-law decay of the radiation. 
Leaver [Q provided a detailed analytical picture of the 
radiative evolution and introduced much of the mathe- 
matical framework adopted in subsequent work. Gund- 
lach. Price, and PuUin ||] provided the most complete 
picture to date; in addition to describing how the radi- 
ation behaves at a fixed distance from the black hole, 
they also described the behavior at the event horizon 
and at future null infinity. Furthermore, they numeri- 
cally simulated the collapse of a self-gravitating scalar 
field, and showed that the power-law tails are a generic 
feature of radiative decay, whether a black hole forms or 
not. These numerical results were reproduced by Burko 
and Ori M . Additional analytical insights can be found 
in Refs. pfli E^nd various generalizations of the basic 
scenario have been considered in pO| , pl| . 

The works reviewed above were entirely concerned 
with gravitational collapse to a non-rotating black hole. 
It is only recently that the analysis was extended to the 
case of rotating black holes. The first authors to con- 
sider radiative decay in Kerr spacetime were Krivan, La- 
guna, and Papadopoulos (l2| (see also [|l3[)i who showed 
numerically that the radiation falls off according to an 
inverse power law. This inverse power-law behavior was 



confirmed analytically by Hod ||Ij] , as well as Barak and 
OriH 

At this point, we have a clear physical picture of how a 
radiative field behaves during gravitational collapse to a 
black hole; the picture is essentially the same for scalar, 
electromagnetic, and gravitational radiation, and it is 
valid for rotating and non-rotating black holes. There 
are three stages to the field's dynamical evolution. At 
the onset of collapse {t = 0), an outburst of radiation 
is emitted; most of this radiation propagates (with dis- 
tortions) directly to infinity. After this initial outburst, 
the field oscillates with frequencies and damping times 
characteristic of the central black hole. This behavior is 
associated with the field's quasi-normal modes of oscil- 
lation, which decay over a time comparable to 100 r^/c, 
where r^ is the event-horizon radius and c the speed of 
light. Finally, these oscillations give way to the inverse 
power-law behavior described above. If the radiative field 
<& is observed at a fixed radius r and the field is static 
prior to collapse, then $ ^ ^-(2/+2) ^g ^ ^ qo, where I is 
the multipole order. 

Analytical |H,R-pl and numerical Q studies of radiative 
dynamics in black-hole spacetimes have revealed that the 
inverse power-law behavior is not sensitive to the pres- 
ence of an event horizon. In fact, power-law tails are 
a weak-curvature phenomenon, and it is the asymptotic 
structure of the spacetime at radii r 3> Tg which dictates 
how the field behaves at times t ^ rg/c. This observation 
begs the question: how is the field's evolution affected if 
the conditions at infinity are altered? For example, what 
happens when the black hole is immersed in an expand- 
ing universe |16[| ? 

This question was first addressed by Brady, Chambers, 
Krivan, and Laguna (hereafter BCKL) in Ref. ||l^. They 
considered the dynamical evolution (both linear and non- 
linear) of a scalar field in Reissner-Norsdtrom-de Sitter 
spacetime, which represents a spherical, charged black 
hole immersed in de Sitter space — an exponentially ex- 
panding universe with a cosmological horizon at a fixed 
radius re- BCKL found that the scalar field decays 
not as an inverse power of time, but exponentially at 
times t ^ rc/c. Their numerical results were compatible 
with the formula $ ~ ^-i^at^ where the decay constant 
Kc ~ c/rc is the surface gravity of the cosmological hori- 
zon. (The surface gravity is precisely defined in Sec. II.) 
For ^ = 0, BCKL found that the field does not decay, but 
that it settles down to a constant, nonzero value. There- 
fore, the usual power-law scenario does not necessarily 
survive when the conditions at infinity are altered. 

Our purpose in this paper is to extend the initial work 
of BCKL; we provide additional details, generalize their 
discussion, and offer new analytical insights. 

We study the dynamics of a scalar field <I> in 
Schwarzschild-de Sitter (SdS) spacetime, which describes 
an electrically neutral black hole (of event-horizon ra- 
dius Te) immersed in an exponentially expanding uni- 
verse (with cosmological- horizon radius Vc). The scalar 
field satisfies the wave equation 



(n-ei?)$ = 0. 



(1.1) 



where □ is the curved spacetime d'Alembertian operator, 
R the spacetime's Ricci scalar, and ^ a tunable, nonneg- 
ative coupling constant. This equation is more general 
than the one considered by BCKL, who limited them- 
selves to minimal coupling, i.e., ^ = 0. We shall see that 
adding this dimension to the parameter space greatly en- 
hances the number of possible late-time behaviors for the 
scalar field. Additional details regarding the SdS space- 
time and the wave equation ( |l.lD are presented in Sec. II, 
which also des crib es the numerical methods employed to 
integrate Eq. (1.1). 

The analysis presented in this paper establishes the 
following physical picture. The SdS spacetime comes 
with two distinct time scales, Tg/c and Tc/c; for the pur- 
pose of this discussion we assume that r^ <C re so that 
the time scales are cleanly separated. The picture elab- 
orated previously for black holes in asymptotically-fiat 
spacetimes — an initial outburst followed by quasi-normal 
oscillations followed by power-law tails — continues to 
hold at times t ^ Tc/c. In fact, over this time scale, 
the field propagates as if it were in asymptotically-fiat 
Schwarzschild spacetime. Conditions at radii r c:i Vc have 
no infiuence on the field's evolution. At later times, how- 
ever, when t becomes comparable to rc/c, the asymptotic 
de Sitter structure starts to play a role, and the inverse 
power-law behavior gives way to a faster, exponential de- 
cay. At times t ^ rc/c the conditions at radii r ^ re no 
longer have any influence, and the scalar field goes as 
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(1.2) 



Equation (L2) reduces to p = ^ when ^ = 0, in agreement 
with the BCKL results. In this expression, notice that 
the square root becomes imaginary when £, > £,c = jq- At 
^ = ^c, therefore, the scalar field undergoes a transition 
at which the pure exponential decay observed for S. < S.c 
becomes oscillatory. This surprising property of wave 
propagation in SdS spacetime was first revealed in our 
numerical simulations (presented in Sec. Ill), and then 
established analytically (as presented in Sec. IV). 

Our final section (Sec. V) is devoted to a discussion of 
the late-time decay of electromagnetic and gravitational 
radiation in SdS spacetime. If $ designates either one of 
these radiative fields, we argue that <I> ~ e~^'*'=* at late 
times, where p is now given by 



p=l+l+0 



(1.3) 



Thus, both electromagnetic and gravitational radiation 
behave, at late times, as a conformally invariant (^ = i) 
scalar field. 



II. EQUATIONS AND NUMERICAL METHODS 

The metric of the Schwarzschild-de Sitter (SdS) space- 
time is given by 

ds^ = -fdt^ + f-^dr^+r'^{d9^+siTi^ed(t>'^), (2.1) 

where 



2M r^ 

J ^ 97 



(2.2) 



with M denoting the black-hole mass; a^ is given in terms 
of the cosmological constant A by a^ = 3/A. We use 
geometrized units such that G — c — 1. The spacetime 
possesses two horizons: the black-hole horizon is at r = r^ 
and the cosmological horizon is at r = re, where re > r^. 
The function / has zeroes at r^, r^ and rg = —{re + re). 
In terms of these quantities, / can be expressed as 



1 



/ = ^- (j' - re){rc - r){r - tq). 



(2.3) 



It is useful to regard r^ and re as the two fundamental pa- 
rameters of the SdS spacetime, and to express M and a^ 
as functions of these variables. The appropriate relations 
are 



2 2 I 1 2 

a = re + r^rc + re 



and 



2Ma'^ = rere{r 



(2.4) 



(2.5) 



We also introduce the surface gravity Ki associated with 
the horizon r — ri, as defined by the relation Ki — 
^\df / dr\r=ri- Explicitly, we have 

(2.6) 

(2.7) 

«o = - " 2';' ^ " ■ (2-8) 

2a^(-ro) 

These quantities allow us to write 
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(2.9) 



and to express the transformation between r and the 
"tortoise coordinate" r* = J f^^ dr as 



r* = ln( — -1) ln(l- — 



1 , / ^ , 
H In 1 

2ko V^o 



(2.10) 



In terms of t and r* we define the null coordinates u — 
t — r* (retarded time) and v = t + r* (advanced time) so 



that the (future) black-hole horizon is located at w = 00, 
and the (future) cosmological horizon is at w = 00. 

We consider a massless scalar field $ in the SdS space- 
time, obeying the wave equation 



(D - f i?)<i> = 0, 



(2.11) 



where □ = .g"^VaV;3 is the d'Alembertian operator, 
R = 12/a^ the Ricci scalar, and £, a curvature-coupling 
constant that we take to be nonnegative. If we decom- 
pose the scalar field according to 



$ 



J2-^i{t,r)Yi„ 



then each wave function ■0; satisfies the equation 

at2 



-in^+Q^]Mt,r)^mr)Mt,r), 



where the potential function is given by 

'1(1 + 1) 2M 2(6^-1) 



Vi{r)^f 
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(2.12) 



(2.13) 



(2.14) 



Using the null coordinates u and u, Eq. ( 2.13| ) can be 
recast as 
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ipi{u,v) = Vi{r)^i{u,v), 



(2.15) 



in which r is determined by inverting the relation r* (r) = 
i(u — u). Our numerical methods are based on this form 
for the reduced wave equation. 



In order to find a unique solution to Eq. (2.15), initial 
data must be specified on the two null surfaces u = 
(say) and v = (say). Because the late-time behavior of 
the wave function is largely independent of the choice of 
initial data, we set tpii^, u = 0) = and use a Gaussian 
profile, centered on Vc and having width cr, on m = 0: 



ipi{u = 0, w) = exp 



{V - Vef 



2a' 



(2.16) 



We numerically integrate Eq. ( 2.15| ) using the finite- 
differencing scheme suggested by Gundlach, Price, and 
PuUin [pj, in which the coordinates u and v increase by 
discrete units A. In the discrete space, the differential 
equation becomes 

i,i{N) = ^i{W) + ^i{E) ~ ^{S) 

- ^ViiRe) \MW) + ME)] + 0(A*), (2.17) 

where we have defined the points N : (u + A,v + A), W : 
{u + A,v), E : [u,v + A), and S : {u, v). The potential is 
evaluated at the central radius Re corresponding to the 
off-grid point {u + ^A,v+ ^A); thus, r*{Re) = ^{v~u). 
The computationally expensive part of the calculation 
resides in the inversion of the relation r*(r), which is 



required to evaluate the potential function Vi{r). For 
values of r* such that 2Ker* < —1, we numerically solve 
the equation 



r 



1 + e^^'"'' exp 
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(2.18) 



by iterations. When 2ncr* > 1, we use 



— = 1 - e" 



exp 



— In 1 +— In 1 



(2.19) 

instead. For all other values of r* , we invert the rela- 
tion r* (r) using the Newton- Raphson method, as imple- 
mented in the Numerical Recipes routine rtsaf e |l8[ . 

Once the integration is completed, the values 
^i{um!i.^,v) are extracted, where w,nax is the maximum 
value of u on the numerical grid. If M,nax is suffi- 
ciently large, then these values give a good approxima- 
tion for the wave function at the event horizon. Similarly, 
for Wmax sufficiently large, the values ipiiu^Vmax) give a 
good approximation for the wave function on the cos- 
mological horizon. Finally, the values of "0; on the line 
u ^ V — 2k are extracted, and expressed as a function of 
t — ^{u + v) — V — k; here. A: is a constant representing 
the value of r* at which the field is evaluated. 



III. NUMERICAL RESULTS 
A. Behavior at early and intermediate times 

Our first goal in this section is to establish that at 
times i <C Tc, the scalar field behaves as if it were in 
Schwarzschild spacetime, and that at times t ^^ Vc, a 
transition occurs in which the inverse power-law decay 
gives way to an exponential decay. The late times t^Vc 
are considered in the next subsection. 

Figure 1 displays the behavior of the wave function 
ipi{t,r) at a fixed radius r, for early and intermediate 
times. The plots show how the pure Schwarzschild be- 
havior (initial outburst, quasi-normal oscillations, and 
power-law decay) is distorted by the de Sitter structure 
of spacetime at large radii. We consider a SdS space- 
time with parameters r^ = 1 and re = 2000, so th at the 
time scales are well separated. We integrate Eq. ( 2.15 ) 
for I = and 1 = 1, setting ^ = 0. We compare the 
field's evolution in SdS spacetime to what it would be in 
pure Schwarzschild spacetime (also w ith r e = 1). In both 
cases the initial data is given by Eq. ( ^.16 ), with Vc = 10, 
(7 = 3, and the wave function is evaluated at r* = 10. 







Schwarzschild-de Sitter (UO) 

Schwarzschild (UO) 

Schwarzschild-de Sitter (kl) 






0.0001 


-*ll ^^"'-^. i 








1e-06 


"''--, 
■-'., 


/^--.. 


1e-08 


"-•>-, 

'-.. 




1e-10 


"-.,_ 

'"""'^--. 


1e-12 





FIG. 1. Absolute value of the wave function i/);(i, r) as a 
function of time t, evaluated at r* — 10 in Schwarzschild 
spacetime (re = 1) and SdS spacetime (re = 1 and re = 2000). 
The cases / = and / — 1 are considered, and the wave 
functions are plotted on a log-log scale. In such a plot, a 
straight line indicates power-law behavior, and a change of 
sign in the wave function is represented by a deep trough. 
We see that the early portion of -i/;! is oscillatory, and that for 
SdS spacetime, tpo changes sign at i ~ 260. 

The plots show that the physical picture given in Sec. I 
is accurate. At early times, the Schwarzschild and SdS 
wave functions display identical behavior, and deviations 
become apparent only when t becomes comparable to Tc. 
For / = 0, the Schwarzschild behavior ipo ^^ t~^ gives 
way to the wave function changing sign at t ^ 260, and 
settling down to a constant value (equal to —1.1 x 10~^) 
at late times. For / = 1, the Schwarzschild behavior 
ipi ~ t^^ is replaced by a faster decay, which eventually 
becomes exponential. Examination of the field's behavior 
at the black- hole and cosmological horizons, where ipi is 
a function of u and v, respectively, reveals precisely the 
same physical picture. 



B. Behavior at late times 

We now consider the behavior of the scalar field at 
times t ^ re- To ease access to these late times, we 
consider a SdS spacetime with parameters Te = 1 and 
Tc = 100, although this no longer offers a clean separation 
of the time scales. We begin by integrating Eq. ( 2.15| ) 
with ^ = 0, for the cases I = 0,1, 2. Once again we set 
Vc = 10, (7 = 3, and we evaluate the wave function at 
r* = 10. 

Figure 2 displays the behavior of the wave function 
ipi{t,r) at late times. The plots reveal that for I = 0, 
the wave function eventually settles down to the constant 
value — 4.4x10"'^. For I y^ 0, the field decays to zero, and 
our numerical results are compatible with the formula 



Ut) 



-iKrJ 



(3.1) 



first discovered by BCKL. The numerically-determined 
values for the decay constants agree with Ikc to within 
1%. Multipoles corresponding to ^ > 3 could not be 
examined because the field decays too rapidly and the 
numerical integration quickly becomes noisy. The late- 
time behavior of the wave function at the blac k-ho le and 
cosmological horizons is also described by Eq. (3.1), pro- 
vided that t is replaced by v on the event horizon, and 
by u on the cosmological horizon. 




FIG. 2. Absolute value of the wave function tpi{t,r) as a 
function of time t, evaluated at r* = 10 in SdS spacetime 
{re — 1 and re — 100). The cases Z = 0,1,2 are considered, 
and the wave functions are plotted on a semi-log scale. In 
such a plot, a straight line indicates exponential behavior. 
Notice that the final change of sign of the wave function occurs 
at t ~ 50 for ; = 0, t ~ 190 for Z = 1, and i ~ 220 for 
I = 2. Notice also that the numerical integration becomes 
noisy when \ipi\ drops below 10"^'*. 

A rich spectrum of late-time behaviors is revealed 
when ^, the curvature- coupling constant, is allowed to 
be nonzero. This dimension of the parameter space was 
not explored by BCKL. Figure 3 displays the late-time 
behavior of ipo for several values of ^, using the same pa- 
rameters as before. For ^ smaller than a critical value ^c, 
the field decays monotonically with a decay constant that 
increases with increasing ^. When ^ > ^c, however, the 
wave function oscillates with a decaying amplitude. As £_ 
is increased away from the critical value ^c, the frequency 
of the oscillations increases, but the decay constant stays 
the same. Similar behavior was observed for I ~ 1 and 
1 = 2, and our numerical results are compatible with the 
formula 



V' 



-ll+g(0]Kat 



where the function g(^) increases with ^ for ^ < ^c, and 
is complex for ^ > ^c- From our numerical results, we are 
able to determine that 0.18 ;$ Cc ^ 0.19. The analytical 



calculations presented in Sec. IV show that the function 
g{S,) is in fact given by 



5(C) 
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1^9^481 -HO(^), 



(3.3) 



so that ^c = js = 0.1875. This relation explains the 
obser ved featur es, a nd our numerical results agre e with 
Eqs. (|U) and (|J) to within 1%. Equation (|]|), with 
t replaced by either v ot u, describes also the behavior 
of the wave function on the black-hole or cosmological 
horizons, respectively. 
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FIG. 3. Absolute value of the wave function i/)o(i, r) as a 
function of time t, evaluated at r* = 10 in SdS spacetime 
{re = 1 and re = 100). Several values of ^ are considered, in 
the interval between ^ = and C = 5- The wave functions 
are plotted on a semi-log scale. The noteworthy features are 
these: (i) For ^ < ^c, the wave function decays exponentially, 
with a decay constant that increases with increasing 5; (ii) for 
^ > ^c, the wave function still decays exponentially, but with 
a decay constant that no longer varies with ^; (iii) for ^ > ^c, 
the wave function oscillates, with a frequency that increases 
with ^. 



According to the physical picture presented in Sec. I, 
the field's evolution should proceed as if the field were 
in pure de Sitter spacetime when t ^ r^- In particular, 
Eqs. (3.2) and (3.3) should describe the late-time behav- 
ior of the wave function in de Sitter spacetime. To test 
this hypothesis, we performed a numerical integration of 
Eq. ( 2.15 ) in pure de Sitter spacetime |1^ with a cos- 
mological horizon at Tc = 1, and initial data of the form 
( ^.16 ) with Vc — I and cr = 0.1. The results for Z = are 
displayed in Figure 4. All the late-time features seen in 
SdS spacetime are reproduced, and we have verified that 
t he field's evolution is well described by Eqs. (3.2) and 
(3.2) (Q, with Te = and k^ = l/^c 



1, to well within 1%. 
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FIG. 4. Absolute value of the wave function tpo{t,r) as a 
function of time t, evaluated at r* = 0.5 in pure de Sitter 
spacetime (with parameter re — 1). Several values of ^ are 
considered, in the interval between ^ = and C = |- The 
wave functions are plotted on a semi-log scale. These plots 
show the same features as in the preceding figure. 



IV. ANALYTICAL RESULTS 



Vir) = f 



'1(1 + 1) 



2(6^-1) 



(4.2) 



We suppress the Z-label to simplify the notation. In 
de Sitter spacetime, the relation between r and /'* = 
J f^^ dr is given by 



„2r* 



1 



(4.3) 



We wish to integrate Eq. ( 2.13 ) starting with a suitable 
set of initial conditions at f == 0. It follows from a simple 
application of Green's identity that at times t > 0, the 
wave function can be expressed as 



i^it.r) 



i:{0,r')g{t;r,r')+^{0,r')g{t;r,r'] 



dr' 
T' 

(4.4) 



Here, ip{0, r') and -0(0, r') are the initial data (an overdot 
indicates time differentiation p^), and g{t;r,r') is the 
retarded Green's function for the reduced wave equation 
( p.l3D ; we also use the notation /' = 1 — r'^. T he G reen's 
function is zero for i < 0, and it satisfies Eq. (2.13) with 
a term fS{t)S{r — r') added to the right-hand side. 



The observation that the late-time behavior of the 
scalar field is dictated by the de Sitter structure of space- 
time at la rge radii is the key to the derivation of Eqs. ( |3.2| ) 
and (3^). Motiva ted b y this observation, we integrate 
the wave equation (2.13) in pure de Sitter spacetime, and 
show that the late-tinre behavior of the wave function is 
indeed given by Eq. (^, with g(^) = | - i^S - 48^- 
Apart from an approximation of late times {t 3> fc), our 
calculations are exact. That such a treatment is afforded 
relies on removing the black hole from the spacetime. 
While this removal produces a large effect at early times, 
the preceding results indicate that the influence of the 
central mass is negligible at late times. (We will have 
to qualify this statement below.) The consid erab ly more 
difficult job of analytically integrating Eq. ( |3.2| ) in full 
SdS spacetime will be undertaken elsewhere [EOf. 



A. Equations 



We consider a scalar field $ in de Sitter spacetime. 



whose metric is of the form (2.1) with 



/ = 1- 



(4.1) 



where we have set re — a = 1 without loss of general- 
ity; thus, the surface gravity of the cosmological hori- 
zon is Kc — I. The scalar field obeys the wave equation 
(bjT ), an d the reduced wave function ■(/'(t, r) defined by 
Eq. ( 2.12 ) satisfies Eq. ( ^.13| ) with the potential 



B. Fourier transform of the Green's function 

Since the late-time behavior oi'4'{t, r) is entirely deter- 
mined by the lat e-time behavior of g{t; r, r'), the deriva- 
tion of Eqs. (3^) and (3_^) begins with a calculation of 
the Green's function. We first consider its Fourier trans- 
form ^(w; r, r'), defined by 



g{t;r,r') 



1 

2^ 



]{uj;r,r')e~ 



■ dw. 



(4.5) 



For each w, this function of r and r' satisfies the differ- 
ential equation 



^+^ -no 



~giu;;r,r') = f5{r-r'), (4.6) 



dr 



where V{r) is given by Eq. ( |4.2| ). It is clear that the 
solution can be expressed in the factorized form 



Kc^;r,r') = -j^0"-'=^(r<)V'"P(r>), 



(4.7) 



where '0'^°s(r) and ■!/'"P(r) are two linearly ind ependent 
solutions to the homogeneous form of Eq. (4.6), with W 
denoting their conserved Wronskian: 



w = ij'' 



dr* dr* 



(4.8) 



We use the notation r< = min(r, r'), r> = max(r, r'), 
and for simplicity, we do not explicitly indicate the de- 
pendence of tp'^'^^, i^^^, and W on uj. 



Equation ( [4.7| ) correctly represents the retarded 
Green's function if '0"P(r) describes waves which are 
purely outgoing at the cosmological horizon; we choose 
the normalization 



r 



1. 



(4.9) 



Also, we must choose for ip^°^{r) a solution which enforces 
the boundary condition ■0'°®(O) = 0; this ensures that the 
scalar field <& is everywhere nonsingular. The calculation 
oig{uj; r, r') therefore reduces to solving the homogeneous 
form of Eq. (4.6) for the two functions ip'^^^{r) and i/j^P{r). 



C. Solutions to the homogeneous equation 

The mathematical problem of finding the solutions 
it^Hji) ^^'^ ip^P{r) has been addressed in the literature 
p2| , P3| . We shall briefly sketch the method of solution. 

The change of variable z = r^ and the factorization 
^ = z"(l - zfX, with a = ^{1 + 1) and /3 = ±|w, 
transforms the homogeneous form of Eq. ( [4.6[ ) to a hy- 
pergeometric equation for X . The fundamental solutions 
are X = F{a±,b±;l + | ; z) , where a± = ^{l + g+ ±iuj), 
b± — ^{l + g^ ± iuj), and 



g±iO = l±lV^~m- 



(4.10) 



Both choices of sign in front of iuj produce a function 
■0(r) that goes to zero (as r'+^) when r ^ 0. Choosing 
the minus sign and the normalization arbitrarily, we set 

V;'-°S(r) =r'+l(l-r2)-"^/2^(a_,6_;? + |;r2). (4.11) 

It is not difficult to show that 

V;"P(r) = 2*"r'+i(l - r2)-'"/2F(a_, 6_; l-icu;!- r^) 



(4.12) 

also satisfies the homogeneous form of Eq. (O), and has 
the asymptotic behavior indicated in Eq. (4.9). On the 
other hand, the solution 

^down(^) ^ 2-*"r'+i(l - r2)'"/2F(a+, b+; l + iuj;l- r^) 

(4.13) 

satisfies ingoing-wave boundary conditions at the cosmo- 
logical horizon: ■0'^°™(?') ~ e"^'^''* when r ^ 1. 

These three solutions are not linearly independent. Us- 
ing Eq. (15.3.6) of Abramowitz and Stegun [|^, it is easy 
to show that 

2— r(/ + |)r(z^) 

V- '^(r) = 7^7 TTTT^ ^ v) 



r(a+)r(6 



r(a_)r(6_) 



V'''°™(r). (4.14) 



This relation, together with the fact that the Wronskian 
between ip'^°'"'^{r) and ?/'^P(r) is equal to 2iuj, imply that 



1 

W 



1 r[i(z + g+-z^)]r[i(r 



9- 



2l+ii^ 



r{l + l)T{l-zu;) 



(4.15) 



D. Pole structure 



If for the moment we ignore the factor l/r(l — iuj), 
Eq. ( 4.15D shows that 1/W [or equivalently g(a;;r, r')] is 
analytic in the complex lu plane, except for isolated poles 
in the lower half plane, at 



lU! 



l + g± + 2n, 



(4.16) 



where g±{£,) is defined in Eq. ( I.ICJ ) and n = 0, 1, 2, • • •. 
There are two poles for each value of the integer n. For 
^ < S,c = Jq, g± are real functions of ^, and the poles are 
located along the negative imaginary axis of the complex 
uj plane. As ^ ^ ^c, the single poles belonging to the 
same n merge, becoming a double pole when ^ = ^c- For 
^ > ^c, the functions g± become complex, and the poles 
move away from the negative imaginary axis. 

The pole nearest to the real axis gives the dominant 
contribution to the Green's function at late times. For 



£, < S,c, this pole is located at iui 

^c, there are two nearest poles at iuj 
For values of £, such that V9 



/ + I-5V9-48C. For 



the double pole is located at iui — I + ^. For £_ > 

/-)-f±W48e-9. 



48^ is an odd integer, 
the numbers u;^ are also integers. This occurs when ^ = 
0, for which the poles are at ito = l+2n and ito = Z-|-3-|-2n. 
It also occurs when £, = h (the value of ^ which makes 
the wave equation conformally invariant), for which the 
poles are at zw = Z -I- 1 -I- 2n and iuj = I + 2 + 2n. These 
are exceptional cases, because most of these poles are 
canceled out by the previously ignored factor l/r{l — iuj). 
This F-function also has integer poles, at iuj — 1 + n' , 
where n' = 0, 1, 2, • • •. It is easy to check that the zeros 
of l/r(l — iuj) cancel out the poles listed previously, so 
that 1/W [or equivalently, g{uj;r,r')] is analytic at these 
frequencies. The only exception occurs for I = when 
^ = 0. In this case, the single pole at iuj — survives. 
This is why the scalar field settles down to a nonzero 
constant at late times. 

The pole cancellation which occurs at integer values 
of iu! has interesting consequences. For example, con- 
sider the case 1 = 1, and values of ^ which are small but 
nonzero. For these values, the pole nearest to the real 
axis is at iuj = 1 + 4^ -|- 0(^^), and this pole exists as 
long as ^ 7^ 0. In the ^ ^ limit, however, the would- 
be pole at iuj = 1 is canceled by the compensating zero 
in l/r(l — iuj). We should therefore expect the Green's 
function — and the scalar field — to display a qualitative 
change of behavior at late times, as the limit ^ = is 



taken. Figure 5 shows that such is indeed the case. It is 
interesting to note that such quahtative changes of be- 
havior do not occur in full SdS spacetime. We return to 
this point in Sec. IV G. 




FIG. 5. Absolute value of the wave function ■0i(t,r) as a 
function of time t, evaluated at r* = 0.5 in pure de Sitter 
spacetime (with parameter Vc — 1). Several small values of ^ 
are considered, together with the special case S, = 0. The plots 
make it quite clear that the late-time behavior associated with 
^ ^ 1 is qualitatively different from the behavior associated 
with ^ = 0. As explained in the text, this qualitative change 
of behavior is caused by pole cancellation in 1/W . 

The presence of poles in ^(w; r, r') makes it very easy to 
calculate the late-time behavior of the Green's function. 
Such a calculation is much more difficult when there are 
no poles, and therefore, this will not be attempted for 
the special cases ^ = (except for I = 0) and £, = ^- In 
the generic situation (all other cases), the contribution 
to the Green's function coming from the pole nearest to 
the real axis is given by 



l{uj;r,r') 



F(iV9-48C) h{r)h{r') 



r{l + ^)T{l-l-g_) lito-LUo)'- 



where 



and 



Hr) 



= r'+i 



(l_7.2)-('+9-)/2 



tLUo^l + ---^9^A8^. 



(4.17) 



(4.18) 



(4.19) 



This calculation involves substituting the asymptotic re- 
lation F(e) ~ 1/e in place of the relevant F-function, and 



evaluating every oth er fa ct or com p rising g{uj 
frequency wq; Eqs. (Il7| ). 
used along the way |25[ . 
(4.19) reduce to 



at the 



( [4. Ill ), ( |4ll , and ( Qsl) a re 
When ^ = Z = 0, Eqs. (tl.lTj)- 



}{uj;r,r') 



rr 

iu> 



(4.20) 



E. Green's function at late times 

The Green's function g(t;r,r') can be obtained from 
g{u;]r,r') b y in verting the Fourier transform; this gives 
rise to Eq. ( |4.5| ). The contour of integration is along the 
real axis of the complex w plane. However, because we 
are interested in the behavior of the Green's function at 
late times, we can safely close the contour in the lower 
half plane without changing the result of the integral. 
The residue theorem then guarantees that g{t,r,r') can 
be expressed as a sum over residues, each pole contribut- 
ing a term proportional to exp{—iuj^t) to the Green's 
function. The dominant contribution comes from the 
pole nearest to the real axis, and we obtain 



g{t;r,r') 



r(lV9^48?) 
F(/ + |)F(l-;-g_) 



h{r)h{r')e 



'\„-ii+a-)t 



(4.21) 



I — ^V9~^48f. Near the cosmological 



where g-{S,) 

horizon, 1 — r^ ~ 4e^^'' , and Eq. ( [4.21|) becomes 



g{t;r~^l,r')^ _2-('+9-) 



F(lV9^48g) 
r{l + l)T{l-l-g^) 



X h{r ) e 



/^p-('+^-)" 



(4.22) 



where u = t — r* . For ^ = Z = 0, both expressions reduce 
to 



l{t;r,r') 



(4.23) 



Equations (4.21)-( 4.23| ) are valid at times i ^ 1, or u ^ 



1. 

We have found that in the generic situation, the late- 
time behavior of the Green's function is given by e~^'^=* 
or e^P'^"'^, where 



1 



p = ; + -- -79^16^. 



(4.24) 



Notice that we have re-introduced the scale parameter 
Kc — ^/rc- For the special case ^ = / = 0, the decay 
constant p is zero, and the Green's function settle s dow n 
to the time-i ndep endent expression given by Eq. ( 4.23 ). 
Equation (4.24) gives the correct late-time behavior of 
the Green's function for ^ < ^c and ^ > ^c- For ^ = ^c, 
the Green's function is obtained by integrating around 
a double pole a.t iuj — I + ^. This requires integration 
by parts, and the end result is that for fixed r < r^ the 
Green's function behaves as {Kct)exp[—{1 + |)Kci]; with 
u replacing t, this is also the correct behavior on the 
cosmological horizon 



By virtue of Eq. (4^), we conclude that in pure de Sit- 
ter spacetime, the late-time behavior of the wave function 
ip{t, r) is also e~^'''=* or e"*"^":", except for the exceptional 
cases — ^ = and £, — ^ — discussed previously. For the 
special case ^ = / = 0, the wave function settles down to 
a constant value. This conclusion is in full accord with 
the numerical results presented in Figs. 4 and 5. 



F. The special case ^ — I — 



G. Analytical results for SdS spacetime 



When C = ^ = 
final value ip (oo 
Eqs. (Q and ( |4.23| ) 
we obtain 



: 0, the wave function settles down to a 
r) which can readily be calculated using 



Re-introducing the scale factor 



-0(oo,r) = -^ 



■ ■ dr' 



(4.25) 



This equation implies that the value of the scalar field 
$ at late times does not vary with r, a necessary conse- 
quence of the cosmological no- hair theorems 26-/ 
also implies that il>{oo,r) scales as 1/rc^ = A/37a prop- 
erty that wa s disc ussed in BCKL. 



It 



Equation (4.25) relates 'ip{oo, r) to the initial data spec- 
ified on the spacelike surface < = 0. We can extract more 
information from this equation, and also test it against 
our numerical results, if we instead express ip{oo,r) in 
terms of initial data specified on the null surface u — 0. 
(We again assume zero data on w = 0.) To effect this 
translation, we assume that the initial data has support 
solely in the "weak-field" region of de Sitter spacetime, 
in which {r/rcY <S^ 1, so that / — 1 and r* ~ r. In 
this region, and in a neighborhood of the surface t — Q, 
we can approximate '(/'(i, r) by the flat-spacetime solu- 
tion ip{t,r) ~ H{v), where H is an arbitrary function 
oi V — t + r. For this solution, -0(0, r) ~ H'{r) (with 
a prime indicating differentiation with respect to the ar- 
gument), and the integral in Eq. (4.25) becomes J rdH. 
After integration by parts, assuming that H(v) has com- 
pact support, we obtain — J H{r) dr, and 



'0(oo, r) 



H{v)dv. 



(4.26) 







Recall that H{v) is the initial wave profile on the surface 
u — 0. For the specific choice (2.16|), we obtain 



ip{oo, r) 



2tt<t 



(4.27) 



We have tested this predictio n aga inst our numerical re- 
sults. The accuracy of Eq. ( 4.27| ) depends strongly on 
the value of Vc, the central position of the Gaussian wave 
packet. If Tc — 1, we find that for Vc < 0.2, our expression 
is accurate to within 1% over a wide selection of values 
for r and a. For Vc > 0.2, the accuracy gets increasingly 
worse as Vc increases; for example, the error grow s to 
15% at Vc — 0.8. In view of the fact that Eq. (4.27) was 
derived under the assumption that the initial data has 
support only in the region r <C ^c of de Sitter spacetime, 
this is the expected result. 

The rate at which the field tp settles down to its fi- 
nal constant value was also explored by BCKL. Their 



numerical results indicated that ijj — ip(oo, r) ~ e 



-QKct 

where q ~ 2. Inspection of Eqs. ( 4.10| ) and ( 4.16 ) reveals 
that the leading contribution to ^ — "0(00, r) comes from 
the 71 = 1 pole in the Green's function, which produces 
q = 2 + g- — 2. Thus, our analytical results are in full 
agreement with the BCKL numerical results. 



Finding the Green's function for a scalar field in SdS 
spacetime is considerably more difficult than what was 
accomplished in this section, and this calculation will be 
the subject of a separate publication |£^. The conclusion 
is that for SdS spacetime, the late-time behavior of the 
wave function is still exponential, tp ~ e"'"^''*, with the 
decay constant p given by the same expression as before, 
apart from corrections of the order of re/rc'. 



p^/+3 ly9^48^+of^ 
2 2 \rr 



(4.28) 



This is the result that was quoted in Eqs. (p]2) and (p.3|). 
In practice, the correction term of order r^/rc is too 
small to be revealed numerically. Nevertheless, its pres- 
ence is important, because it prevents p from ever becom- 
ing an integer. This means that the cancellation of poles 
at integer values of ito, discussed in Sec. IV D, is a phe- 
nomenon that occurs only for pure de Sitter spacetime; 
pole cancellation does not occur in full SdS spacetime. 
Consequently, the cases ■f = and C = ^ are not excep- 
tional in SdS spacetime, and such qualitative changes of 
behavior as depicted in Fig. 5 do not occur. 



We have numerically tested the validity of Eq. (4.27) 
for SdS spacetime. As expected, the accuracy of this for- 
mula is worse than for pure de Sitter spacetime: With 
Te = 1, Tc = 500, a = 10, and evaluating the wave func- 
tion at r* ~ 100, we find that the error is less than 5% if 
100 < Vc < 200, but that it chmbs to approximately 15% 
when Vc ~ 50 or Vr ^ 400. 



V. ELECTROMAGNETIC AND 
GRAVITATIONAL RADIATION 



The scalar- wave equation (1.1) serves as a model for 
the propagation of electromagnetic and gravitational 
waves in the SdS spacetime. In this section we indicate 
how our results can be generalized to cover these physi- 
cally important cases; a more complete discussion will be 
found inRef. (2^. 

It is well known that electromagnetic and gravitational 
perturbations of the SdS spacetime can be analyzed in 
terms of a master potential $ which is related to the 
components of the perturbing fields pa] . After a decom- 
position in spherical harmonics, the wave function ij jijt, r ) 
is found to obey an equation of the form of Eq. (2.13), 
but with the generalized potential 



Viir)^f 



1(1 + 1) 2(s2-l)M' 



(5.1) 



where s = 1 (s = 2) for electromagnetic (gravitational) 
perturbations. In both cases the multipole index / is 
restricted by / > s. 



The considerations of the preceding sections indicate 
that the late-time behavior of "0; is i nsen sitive to the pres- 
ence of M-dependent terms in Eq. (5.1). Removing these 



terms, we obtain the potential of Eq. (42) with £, = ■^■ 



On these grounds, we expect that the late-time beh avio r 
of $ is exponential, <& ~ e~P'^''*, with p given by Eq. ( |l.2| ) 
with^= i: 



p^l+l+O 



This result is borne out in our numerical simulations. 
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